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O 

, Abstract. We study finite-time blow-up for pseudo-relativistic Hartree- and 

^ ' Hartree-Fock equations, which are model equations for the dynamical evolu- 

bJQ' tion of white dwarfs. In particular, we prove that radially symmetric initial 

^ . configurations with negative energy lead to finite-time blow-up of solutions. 

' Furthermore, we derive a mass concentration estimate for radial blow-up so- 

lutions. Both results are mathematically rigorous and are in accordance with 
Chandrasekhar's physical theory of white dwarfs, stating that stellar configu- 
rations beyond a certain limiting mass lead to "gravitational collapse" of these 
objects. Apart from studying blow-up, we also prove local well-posedness of the 
initial-value problem for the Hartree- and Hartree-Fock equations underlying 
' our analysis, as well as global-in-time existence of solutions with sufficiently 

^ small initial data, corresponding to white dwarfs whose stellar mass is below 

QQ ■ the Chandrasekhar limit. 

o 

O ■ 1- Introduction and Description of the Problem 



> 

o 



> 
X 



This paper is a contribution to the mathematical physics of white dwarfs and 
neutron stars. White dwarfs are dense stars composed of electrons and nuclei, which 
»^ I form a completely ionized plasma. The electrostatic Coulomb forces between these 

■ particles establish local electric neutrality to a high degree. For this reason, these 

^ , forces are screened almost perfectly and make only a very modest contribution to 

the dynamical evolution and the energy of a white dwarf. Local electric neutrality 
implies that the spatial and momentum distributions of nuclei are approximately 
. , equal to those of the electrons. Since the masses of nuclei are much larger than 

' the mass of an electron, the leading contribution to the total kinetic energy of a 

white dwarf comes from the electron gas, while the main contribution to its potential 
energy is due to the gravitational interaction among the nuclei. To simplify matters, 
we consider a single species of nuclei of electric charge Ze and mass mz 3> m, where 
— e is the electric charge of the electron and m its mass. Throughout this paper, we 
use units such that Planck's constant h = 1 and the velocity of light c ~ 1. 

1.1. Heuristic Discussion. Let N denote the number of electrons in a white 
dwarf, and let R be its radius. Then the number of nuclei is N/Z, and the av- 
erage momentum, p, of a nucleus or electron is given by p ~ N^^^ /R. A rough 
estimate for the groundstate energy, E{N), of such a star is thus given by 



2 ^^^2 



(1) iJ(iV)=mm{7Vv^^T^+|mz-^|) ^} 

The constant ^mz is the rest energy of the nuclei and will be subtracted from the 
groundstate energy henceforth; the constant G is Newton's gravitational constant. 
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The function of p within curly brackets on the right side of Q has a minimum 
> — oo only if 

(2) 7V<iV„:=(Gmz)"'/'z3, 



where A^cr is the "Chandrasekhar number. " Thus, a white dwarf of total mass 

A', 
Z ' 



M ~ y-mz larger than the so-called "Chandrasekhar mass' 



(3) Me, := ^mz, 

is energetically unstable and is expected to undergo gravitational collapse; see [5]. 

For N < iVcr, the momentum p minimizing the right side of is of the order of 
m; (more precisely, p > m{N /NdY^"^ , with p oo, as iV ^ Nd-—, see e.g. ,14). 
The use of relativistic kinematics for the electrons in a calculation of the groundstate 
energy of a white dwarf is therefore mandatory. Furthermore, for typical white 
dwarfs, the ratio between the Schwarzschild radius, 2GM, and the radius R of the 
star, as determined by is of the order of 10^'', so that effects of general relativity 
are unimportant, and gravity can be described by unretarded Newtonian two-body 
forces. 

The outline of a heuristic description of neutron stars is similar, except that 
effects of general relativity become more important. 

1.2. Hartree- and Hartree-Fock Equations. We now propose to describe a 
white dwarf, or a neutron star, quantum-mechnically, but within the approximation 
described above and assuming that the number N of electrons is conserved. This 
leads us to consider the Hamilton operator 

^ I 1 

(4) if(-) = Ev/^^^-« E y-h^r 

k=l l<k<l<N ' ' 

where pk — —i'^x^ and k — Grn^/Z^, acting on the Hilbert space 

(5) T^W = (l2(m3)®c2)'''^. 

Here L^(R'^) is the space of square- integrable one-electron wave functions on physi- 
cal space K'^, is the space of states of the spin of an electron, and "/\N" denotes 
an TV- fold antisymmetric tensor product, in accordance with the fact that electrons 
(and neutrons) are fermions, i. e., they satisfy the Pauli principle. In the following, 
electron spin plays a completely uninteresting role. To simplify our notation, we 
will therefore ignore it. 

Special state vectors in the Hilbert space Ti^^^ are given by Slater determinants. 



(6) A:=ViA...AVa, 

where ipi, . . . ^ipN are N orthonormal one-particle wave functions; i. e., 

(7) {'4>k,i>i) ^ Ski, for all A:,; = 1,..., TV, 

with (•, •) the scalar product on L^(R'^), (remember that we neglect electron spin). 
We note that the (dimensionless) couphng constant k = Gm^/Z^ is tiny, n ~ 
0(10^'^*), while the number N of electrons in a star is huge: N < 0(10^^). One ex- 
pects that, in this regime, the groundstate of the Hamiltonian i/^^-' is well approxi- 
mated by a Slater determinant, Aq; see |1(JII11) . Moreover, the quantum-mechanical 
time evolution, as described by the one-parameter unitary group {exp{—itH^^^ )}t^R, 
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is expected to evolve a Slater determinant, 0i A . . . A ^^v, describing the state of the 
star at time t — 0, to another Slater determinant, 

(8) VlW A...AVAr(t), 

at time t > 0, up to an error term that tends to in the "mean-field limit" k — > 0, 
with N - Oin"^/^). In ©, the one-particle wave functions il>i{t), . . . ,ipN{t) are 
solutions of the N coupled equations 

N N 

(HF) idti^k = V-A + m^yjk - E (n * i^'HV'fc + J2Mr\* {^'^^1)' 

1=1 1=1 

with initial conditions ^/'^(t = 0) — 4'k, k ~ 1, . . . , TV, and time < t < T, where 
< T < CX3 is the maximal time of existence of the solution {V'fclOl^i- (HF), the 
pseudo-differential operator \/~A + m? , which is defined by its symbol ^/p^ + in 
Fourier space, describes the kinetic energy, including the rest energy, of an electron, 
and the symbol * denotes convolution of functions on R^. The last term on the 
right side of (HF) is the so-called "exchange term, " which is a consequence of the 
Pauli principle. When compared to the second term on the right side of (HF), the 
"direct term," it is subleading in the mean-field limit (N — > oo). It is therefore 
often neglected. Then equation (HF) is replaced by the N coupled equations 

N 



(H) idti'k = V-A + m^^k-J2in* 

1=1 1^1 

with i/jkit = 0) = (pk and k ^ 1, . . . , N. 

Equations (HF) are called (dynamical) Hartree-Fock equations, while (H) are 
called Hartree equations. These systems of evolution equations are the main char- 
acters studied in this paper. 

The logics leading from the quantum-mechanical time evolution generated by 
the Hamiltonian i/^^^ to the nonlinear evolution equations (HF) and (H) for N 
orthonormal one-particle wave functions, ipii't)^ ■ ■ ■ j'4'N{t), in the mean-field limit, 
K with N = 0{k~^^'^), has been studied in 

1.3. Hamiltonian Structure and Conserved Quantities. Let H-^/^(M.^) denote 
the inhomogeneous Sobolev space of index 1/2. We define 

(9) = (i/l/2(K3))XJV_ 

This space can be interpreted as an affine "classical phase space" with complex 
coordinates ^ — (V'l, . . . , V'w) and '4' = {^p^, . . . , V-'jv)- The symplectic 2-form, uj, is 
given by 

. N 

(10) = ^EdV'feAd^fe. 

k=l 

We define a Hamilton functional, 'H^^\ on p(^) by setting 

(11) H^f^*, *)=£#(*), 

with ^ = HF or H, respectively, where the energy functionals ^hf and £h are 
defined in Section|2 below. Then eqs. (HP) and (H) turn out to be the Hamiltonian 
equations of motion corresponding to the Hamilton functionals 'H'^y ''^iT'^ 
respectively. Formally, the quantities 

f#(*(i)) and {'^k{t),Mt)), 
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with fc, Z = 1,...,7V, are conserved under the Hamiltonian flow determined by 
eqs. (HF), (H), respectively. These conservation laws play an important role in our 
analysis. 

1.4. Notation. Throughout this text, we make use of inhomogeneous and homoge- 
nous Sobolev spaces of order s, denoted by H'^{U.^) and ij'*(R^), which are equipped 
with norms \\u\\h- = + V^Y^^uWl^ and = \\{-Ay^'^u\\L2, respectively. 

The scalar product on L^(M^) is defined as {u, v) — J^^ uv dx. 

With some abuse of notation, we sometimes identify collections of wave functions, 
* = {V'fclfcLu "^ith ordered tuples = {tpi, . . . ,i{}n)- For solutions of (HF) and 
(H), this procedure is legitimate, thanks to their U (A^)-gauge invariance; see Section 
13 below. 

In what follows, we write X <Y ii X < CY, where C is some universal constant. 
With regard to physical applications of (H) and (HF), we remind the reader that 
we use units such that h — c— 1. 

2. Main Results 

We begin by reviewing some aspects of equations (H) and (HF). First, we re- 
call from Section 1 that both sets of coupled equations exhibit (formally, at least) 
conservation of energy. That is, the Hartree energy, 

(12) £H(*)=f:(^.,y^^T^V^.)-? / / ^^^^dxdy, 
^ 2 7r3 7r3 \x-y\ 

and the Hartree-Fock energy, 
(13) 

«*) = E(^.,V^^T^^.)-f / / ^^(-W(^)-IP.(^.i/)P d.d,, 

^ Jmdm^ \x-y\ 

are conserved for solutions = {"^fel^i of (H) and (HF), respectively. Here and 
in what follows, we make use of the density matrix, 

N 

(14) p^{x,y) ^^'iljk{x)ipf,{y), 

k=l 

and the particle density, 

(15) p^s,{x) ^ p^s,{x,x). 

In addition to conservation of energy, we also have conservation of the particle 
number (proportional to the stellar mass) given by 

(16) AA(*) = / p^{x)dx. 

As we will see below, the conservation of A/'('4') is a special consequence of the U{N)- 
gauge symmetry; i.e., every transformation ipk '—^ '^i=iTki'4'ii with T G U{N), 
yields another solution of (H) and (HF), respectively. 

2.1. Initial- Value Problem. Our first main result states local well-posedness of 
the initial- value problems for (H) and (HF), provided that the set of initial data 
belongs to H^{M.^), for some s > 1/2. 
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Theorem 1 (Local Well-Posedness) . Let (#) denote either (H) or (HF). Sup- 
pose that s > 1/2 and let N > 1 be an integer. Then the initial-value problem for 
(#) is locally well-posed in H'^{M.^). 

By this we mean the following. For every collection of initial data, ^ — {4>k}k=i 
^^^(R^), there exists a unique solution, = {ipk{'t)}k=i ^ H^{M.^), solving (#) 



holds, for all k = \, . . . ,N . Here < T < oo denotes the maximal time of existence, 
and T < oo implies that limt^T- ||V'fc(0llffi/2 — ^ holds, for some k = I, . . . , N . 
In addition, the solution "^{t) depends continuously on and 

and M{^{t)) = Af{^) 

hold for all times < t < T . Moreover, we have that 

for all 1 < k,l < N and all times < t < T . 

Remark. The proof of this theorem proceeds along the hnes of |S] and wiU only 
be sketched in Subsection 13.11 below. We remark that no use of Strichartz-type 
estimates for the propagator, e^'*^^'^+™^ , is made throughout this proof. Note that 
when considering initial data below the energy norm, e. g., belonging to 7?*(M'^), for 
some s < 1/2, one would have to resort to such estimates; see Since we are only 
interested in finite-energy solutions of (H) and (HF), we have no reason to pursue 
this issue here. 

The next theorem shows that sufficiently small initial data lead to global-in-time 
solutions. The smallness condition corresponds to a number of particles below the 
Chandrasckhar number, iVcr, mentioned in Section 1. 

Theorem 2 (Global solutions for N < A^cr)- Every solution of either (H) or 
(HF), given by Theorem^ exists for all times, < t < oo, whenever the corre- 
sponding initial data, 4> — {<j>k}k=i -ff^(M''), form a collection of L^-orthonormal 
functions whose number N — satisfies 



Here k^-c > is a universal constant of order 1. 

Remark. If, in addition, the initial data $ satisfy f#($) < Nm (the rest energy 
of the electrons) then ||'5'(t)||iP does not tend to 0, as i tends to oo, for any p > 2. 
The physical interpretation of this result is that '4'(t) describes the evolution of a 
bound configuration of matter forming a star-like object. For details and a proof of 
a closely analogous result, see |7]. 

2.2. Finite-Time Blow-Up. We now turn our attention to finite-time blow-up for 
the Hartree equation (H), which, by Theorem |21 is only encountered for sufficiently 
large initial data. 

Theorem 3 (Radial blow-up for (H) with negative energy). Let $ = 

{4'k}k=i C^(R^) be a collection of functions with the property that p^{x) — 
SfcLi is radially symmetric. If the Hartree energy is strictly negative, i. e.. 



such that 



V'fc(O) = 0fc and VfeeC°([0,r);iJ^(R3))nCi([0,T);i?^^i(R3)) 




fH(*) < 0, 
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then the solution, ^{t) — {V'fe(^)}fcLi; of ^) with initial data $ blows up within 
finite time. That is, we have that 

lim^ ||V'fc(i)||_Hi/2 = oo, 
for some k = 1, . . . , N and some T < oo. 

Remarks. 1) It is not difficult to construct an initial configuration €> with £u{^) < 
0, as follows: We consider a ball in M.^ of radius R > centered at the origin. We 
then pick an integer iV > and let 4>i, (t>N denote eigenfunctions of the Laplacian 
with Dirichlet boundary conditions at the boundary of the ball corresponding to 
the lowest A'^ eigenvalues and spanning a rotation-invariant subspace of L^(K'^). In 
accordance with the Pauli Principle, we choose $ := {4>k}k=i- The relativistic 
kinetic energy of this configuration is proportional to A^*/^ , while the gravitational 
potential energy is proportional to — iV^. Thus, if N is sufficiently large then £h(*) 
is strictly negative. 

2) We have also found an analogous blow-up result for (HF), but with the addi- 
tional assumption that each function 4)k{x) has to be spherically symmetric. From 
the physical point of view such a hypothesis appears unnaturally strong, so that we 
refrain from formulating this blow-up result for (HF) as a theorem. 

3) The requirement that 0^ e C^(R3) can be relaxed to weaker conditions on 
regularity and spatial decay. For the sake of simplicity of our presentation, we will 
not pursue this issue here. 

4) By invariance of # = {(t>k]k=i under spatial rotations we mean that, for 
every R e SO{i), (j)k{Rx) = TA=iTki{R)<Pi{x) holds, where T{R) £ U{N) is some 
unitary matrix. This implies in particular that the density matrix obeys p^{x, y) = 
pq,{Rx, Ry), for any R e 5*0(3). Moreover, we remark that it is easy to see that the 
corresponding unique solution, \I'(t), of (HF) or (H) is also invariant under spatial 
rotations, for all times, < t < T, provided that '4'(0) = # has this property. 

Our last result shows that, when approaching the time of blow-up, any radial 
blow-up solution of (H) or (HF) exhibits a concentration of particles at the origin, 
whose number is at least of order of the Chandrasekhar number A'cr- 

Theorem 4 (Chandrasekhar mass concentration for radial blow-up). Let 

^(t) = {i^kit)}k=i '^'^ H^/^-valued solution of either ifflj) or JfHF]) that blows 
up at time T > 0. Moreover, suppose that p-q,(^o){x) — SfeLi IV'fc(0, a^)P is radially 
symmetric. Then, for every R > 0, we have that 



liminf / p-^{t){x)dx> (— — ) 



3/2 



'|2:|<-R 

where p<si(t){x) ~ \4'k{t, x)\^ . Here k„ > is the same universal constant as 

in Theorem\^ 

Remark. In view of the physical interpretation of (HF) and (H) discussed in Sec- 
tion 1, it would be of considerable interest to gain more insight into the properties 
of blow-up solutions for these equations and to arrive at a state of affairs compa- 
rable to what is known about blow-up for nonlinear Schrodinger equations (NLS) 
with L^-critical, focusing nonlinearities; see, e. g., the monograph ^ and references 
given there; (see, in particular, jl2lll5j for mass concentration of blow-up solutions 
for NLS). 
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3. Proof of Main Results 

The proofs of Theorems ^3 are extensions of arguments derived in show- 
ing local and global well-posedness, as well as finite-time blow-up for the pseudo- 
relativistic Hartree equation, i.e., the equation 



(17) idt^lj = V-A + m^^P- {\x\-^ * IV'HV', 

where ijj : [0,T) x R'^ ^ C. We therefore only sketch the proofs of Theorems 

in Subsections I3.1H3.3I In contrast, the proof of Theorem 0] is given in detail in 

Subsection 13.41 

3.1. Proof of Theorem^ For definiteness we consider the initial- value problem 
for (HF) in i?*(R'^), and we observe that all arguments apply to (H), with almost 
no change. 

The initial-value problem for (HF) can be written as follows. 



(IVP) 
Here 



i9t* = V-A + m2 * + kF(*), 
*(0) = * G i/"^^, < i < T. 

TS,N / TjsiTa3\\^N 



(18) H'-^ := {H%M.^))' 
is the A^-fold cartesian product of iJ''(R'^), equipped with the norm 

N 

(19) = 

k=l 

With some abuse of notation, we sometimes identify ^ = {ipk}k=i with the vector 
^, and, likewise, ^ with In IIIVP|) the nonlinearity, F = {Fi, . . . , F^), is given 

by 

N N 

(20) (Ffc(*)) = -Y.{\x\-' * \^Pif)^Pk+J2M\^\'' * (^M). 

1=1 1=1 

Lemma 1. Suppose s > 1/2, and let N > 1 be an integer. Then F : i?"'^ — > H'^'^ 
is locally Lipschitz such that 

(21) ||F(*) - F(*)|U.,« < + - 



(22) ||i?(*)||^.« <||*|i^.,„||#|l^.,«, 

for all #, * e iJ"'^, where r = max{s - 1, 1/2}. 

Remark. The proof of Lemma^for the first term in (|20|l follows from a straightfor- 
ward extension of LemmaS], where an analogous result is shown for the nonlin- 
earity J : H^{M.^) — !■ _ff*(IR^) with J(u) = (|a;|~^ * |wp)w, corresponding to Hartree 
nonlinearities and = 1. Also, the proof of the estimates (|21|l . H22() for the second 
term in l|2U|l . i. e., the "exchange term", can be shown in a similar fashion. 

By Lemma ^ local-in-time existence and uniqueness of ^{t) E H^'^, as well 
as continuous dependence on now follow by standard methods for evolution 
equations with locally Lipschitz nonlinearities; see jH| and references given there. In 
addition, estimate and Gronwall's inequality allow us to deduce that, for any 
s > 1/2, 

sup \\4{t)\\Hs.N <C{T,,\\^{0)\\hs,n, sup \\4{t)\\H./2.N), forO<n<r. 

0<t<T. 0<t<T. 
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In particular, this implies that the maximal time of existence of any _ff*-valued 
solution of (HF), with s > 1/2, coincides with its maximal time of existence when 
viewed as an 7J^/^-valued solutions; see also [S]. 

Finally, we note that [/(A^)-charge conservation is due to ^(■(/'^(i), ^/(t)) = 0, 
which follows from a direct calculation. Moreover, conservation of energy stems 
from the fact that ■^£{'^{t)) — holds whenever the initial data satisfy $ C H^'^ . 
To prove conservation of energy for i?''(]R^)-valued solutions when 1/2 < s < 1, one 
can proceed in a standard way, i. e., by using the continuous dependence on initial 
data in H"'^ and by appealing to the density of H^'^ C H'^'^ . This completes our 
sketch of the proof of Theorem ^ ■ 

3.2. Proof of Theorem m Let (#) either stand for (H) or (HF). Then, by our 
hypothesis on the initial data and Theorem ^ we have that the solution, ^{t) = 
{-(/'fc(0}fcLi C if^/^(R^), of (#) forms a collection of L^-orthonormal functions. By 
Lemma 2 (in Appendix A) and the estimate 

(23) / / ^i^dxd2/<C(/ p{x)dxf\l p{xY''dx\ 

for some constant C ~ 1 (see JT] and reference given there) , we obtain that 

N 

ll,/„,.mi|2 _ _ 

2 7r3 7r3 \x~y\ 



£#W>2^\\Mt)\\Hi/2-i^ I r— n dxdy 



fc=i 



N 



(24) >[l-— (/ P^w{x)dxf']J2\\Mt)\\%,., 

L '^cr JR3 

for some universal constant Kcr of order 1. Using that £#('$'(/;)) = £#{^) and 
/k3 = we deduce the a-priori bound 

N 

(25) sup \\Mmlu2 <N+{1- ^N'/^)-'£#{^), 
o^*<'^ k~i 

provided that N < {ncv/^^)'^^^ holds. By Theorem^] this implies that the maximal 
time of existence is T = cxd. ■ 

3.3. Proof of Theoreml^l The proof of TheoremOfoUows 6 , with some cosmetic 
changes only. We remark that smoothness and sufficient spatial decay of the initial 
data (e.g., that (j)k & H^(M.^) and \ {4>k, \x\'^4'k) \ < oo) guarantee that all quantities 
involved in the following calculations are well defined. 
First, we notice that 

N 

(26) a{t) ■.^J2{Mt),Mk{t)), with A:--i(x- V + V-a;), 

fc=i 

is found to satisfy the differential inequality 

(27) d{t) < fH(*(t)) = fH(*), for 0<t<T. 
Further, a calculation similar to the one in (6j shows that 

N 

(28) m{t) := ^iipkit), Mipk{t)), with M := x ■ V-A + m2x, 

k=l 

obeys 

(29) m{t) <2a{t) + Ci, for < t < T, 
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where Ci > is some constant only depending on Af{^). In our proof of H29|l . 
we make use of Newton's theorem for radially symmetric densities p*(t)(a;), which 
forces us to assume that p^{x) be radially symmetric; see also Remark 4) following 
Theorem 13 

By combining 1271) and (|29|l and integrating, we thus obtain 

(30) m{t) < fH(*)t^ + Cit + Ca, for < t < T. 

Since m{t) is a nonnegative quantity, we conclude that if £B_{'^{t)) — £-r{^) < 
then the maximal time of existence T must be finite. By Theorem^ the solution, 
\I'(t), of (H) thus blows up at time T, which is bounded from above by the positive 
root of the right-hand side of H3U|) . ■ 

3.4. Proof of Theorem |31 We consider (HF); but the proof for (H) is almost 
identical to the one presented below. 

The proof makes use of variational arguments; see also |12[ [T5] for a proof of 
mass concentration for NLS. We argue by contradiction, as follows. Suppose the 
conclusion of Theorem^were not true. That is, there exists some R > such that 

(31) lim / p^(f^)(x) da; < ( — 

holds, where {tnjnsN is a sequence of times increasing to T, as n ^ oo. Here 
Kcr > is the same universal constant as in Theorem |5| 

To prove that assumption 1311) cannot hold, we begin by introducing a collection 
of N sequences, ~ {V'fe^nl^^i; where 

(32) iikA^) ■.^<jit„)-^/^Mtn,<j{tn)~^x). 

Here the strictly positive function a{t), defined on [0, T), is given by 

N 

(33) a(t) :-^||Vfc(i)||^,/.. 

k=l 

Note that, since ||'0fc(OllL2 = ||V-'fe(0)||L2, we have that cr(t) oo, as t — > T— . Next, 
we define the functional 



2 ^ f f P*(2;)p*(2/) - |p*(a;,?;)| 



(34) £(*) 110,11^,,, / / dxdy, 

^ « ' 2 7r3 7r3 \x - y\ 

for any collection $ = {'Pk}k=i -ff^^^(K'^). A simple calculation then yields 

(35) f (*„) = cT{tnr^E{^{tn)). 
Moreover, by using conservation of mass and energy, it follows that 

(36) l^(*(in))| < |f (*(0))| + mAA(*(0)) 



where we use that < (0, V— A + m?4>) — ||0||^i/2 — ^'^II'/'IIls- Combining estimate 
()36|l with the fact that a{tn)^^ ^ 0, as n ^ oo, we conclude that 

(37) lim f (*„) = 0. 

n — >oo 

For later use, we record that this implies that 

(38) lim / / p~.S-)p~.Sy)-\p~.S-^y)\\^^^^i^ 

"^°°Jr3Jr3 \x-y\ K 

since Y^^=i llV'fe.nll^i/a = 1 for all n e N, by construction of i)k,n- 

As a next step, we notice that = {V'fc,n}^i is a collection of bounded se- 
quences in i?^/^(R'^). Hence, by passing to a subsequence, we find that ■(/'fc,™ ^ V'fe,*i 
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weakly in H^^^{M.^), as well as ipk,n{x) ^ ''Pk,* (x) pointwise, for almost every x € R'^, 
as n — !■ oo. Correspondingly, we write yP* = {V'fc,*}fcLi: ^^id the sequence of density 
matrices obeys {x,y) Pq, {x,y) pointwise, for almost every x and y in R^, as 
n — > c». Next, we claim that 

(39) /0^„(a;) PiaS^)^ strongly in iP(R3), for 1< p < 3/2, 

after possibly passing to a subsequence. To prove H39(l . we define the functions 



N 



\ k=l 



(40) /t„(x) := ^p^Ja;) = 
One can show (using ^ Thcorm 7.13]) that 

N 

(41) llAJnII?,!/. <5]||V^fc,„|lW- 

fe=l 

Thus, {KnjneN forms a bounded sequence in H^/'^{W?) and, by hypothesis of The- 
orem 0] and Remark 4) following Theorem ^ we actually have that {KnlneN C 
H^l^iM?), which denotes the subspace of spherically symmetric functions in H'^l'^iJS?) 
Since the embedding 

(42) L^i^') 

is compact if and only if 2 < p < 3 (see J3]), we deduce (after passing to a 
subsequence) that k„ Kcr, as n — s- oo, strongly in LP{R^), for every 2 < p < 3. 
This proves our claim H39|l . 

Combining all the convergence properties shown above, we infer that 

f f p^^{x)p^ {y)~\p^ {x,y)\'^ 
hm / / ■ ■ dxdy 

f f p^M)PiSv)-\PiS^^y)\\ J 

(43) = / / ■ ■ dxdy. 

Here we have used H39|l together with the Hardy-Littlewood-Sobolev inequality, as 
well as the dominated convergence theorem combined with the pointwise estimate 
\Pq, (2^1 y)P ^ P^ {x)p^ (y), which follows from the Cauchy-Schwarz inequality. 
We note, further, that 

(44) 0= lim £'(*„)>£(*,), 

n — >oo 

by using H43|l and the weak lower semicontinuity of the first term in (|34|l . In addition, 
we obtain that < (V'fc,* , "0;,*) ^ ^ki (in the sense of Hermitian N x A^-matrices), 
since i/'fe.n ^ V'fc,* weakly in L^(R'^). Invoking Lemma 2 (see Appendix A) for 
= {V'fc,*}^! we find (similar to the proof of Theorem|2) that 

(45) e{^_^)>k\i~—{[ p^{x)dxf'^\[ p^{xf'^Ax, 

where K > 1.63 is the constant from Lemma 2, and Hcr > denotes the same 
universal constant as in Theorem |3 Moreover, we deduce from H43I) and (|38|l that 
p^ (x) ^ must hold. Thus, we see that inequalities (|45(l and l|44|l imply 

(46) / p^ ix)dx>(^f\ 



DYNAMICAL COLLAPSE OF WHITE DWARFS 



11 



By using (|46|l . we find that assumption H31|) leads to a contradiction as follows. 
In view of we deduce that {x) {x) strongly in Lj'Qj,(M^), as n ^ oo. 

For every A > 0, we thus obtain that 



/ (x) da; = hm / (x) dx 

J\x\<A * "^°°,/|a;|<A " 



= lim / p^^t„)ix)dx 

J\x\<a{f,,)-^A 



(47) < liminf / p^u ){x)dx, 

J\x\<R 

where we use that a{tn)^^ — > as n oo. Since estimate (|47|) holds for every 
A > 0, we deduce from H46|l that assumption leads to a contradiction. This 
completes the proof of Theorem 01 ■ 

Appendix A. Lower bound for Kinetic Energy 
The following result is a slight extension of an estimate derived in 0]. 

Lemma 2. Suppose * = {(l)k}k=i C H^^'^{M.^) satisfies < {4>kAi) < ^ki in the 
sense of Hermitean N x N -matrices. Then 



N 



>k( p^{xf'^dx, 

JR3 



where p^{x) — X^a^i l'/'fc(a;)P and K > 1.63 is some constant. 

Proof. First we remark that both sides of the inequality to be shown are invariant 
under unitary transformations of the (^'s, i.e., under the transformations (pk ^ 
TldLi ^ki4>U for a.n arbitrary unitary matrix A £ U{N). Therefore we can assume 
without loss of generality that 

(48) {(f>k,(t)i) = XkSki, with < Afc < 1, 

where we have also discarded any possible zero vector, (j)k = 0, corresponding to 
Afc = 0. 

Next, we consider the N x A^-matrix, H ~ {hki)i<k.i<N , with entries 

(49) hki = (0fe, V^^i) - c{q}k,UPi). 

Here U{x) p$(a;)^/'^, and c > is some constant to be chosen below. Since H is 
Hermitean, there exists B G U{N) such that H = B*HB has entries 

(50) hki = ^kSki, 

with eigenvalues ei < £2 < • • • < ^n- Let {Ej}j>o denote the set of negative 
eigenvalues of the relativistic Schrodinger operator V— A — cU, acting on L^(M^), 
and we consider the set of orthogonal vectors given by {ipk ■ (pk — J2i BkiPh ?fe < 0}. 
Noting that < {(pk, Pk) < 1 holds, we deduce that 

N N N N ~ 

J2hkk = J2^k>J2^k>Yl ( ,t\ \ ) 

fe=l k=l fc=i fc=i \<Pk,<Pk) 

(51) > V > -Lc^ I U{xf dx. 

Here the first inequality in the second line follows from the min-max principle 
applied to the set of orthonormal vectors {^fc/ll^fellL^} and the operator ^—A — cU. 
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Moreover, the last inequality in H51|) is a standard estimate, where L > is some 
constant; see 0|. By choosing c = 2^'^I'^L''^I'^ ^ we complete the proof of Lemma|2 
where the lower bound K = 2i/A- 2^'^I'^L~^I'^ > 1.63 follows from known boimds for 

L. m 
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